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ABSTRACT 



Context. The APEX Telescope Large Area Survey of the Galaxy (ATLASGAL) is the first systematic survey of the inner Galactic 
plane in the sub-millimetre. The observations were carried out with the Large APEX Bolometer Camera (LABOCA), an array of 295 
bolometers observing at 870 /jm (345 GHz). 

Aims. Here we present a first version of the compact source catalogue extracted from this survey. This catalogue provides an unbiased 
database of dusty clumps in the inner Galaxy. 

Methods. The construction of this catalogue was made using the source extraction routine SExtractor. We have cross-associated 
the obtained sources with the IRAS and MSX catalogues, in order to constrain their nature. 

Results. We have detected 6639 compact sources in the range from 330 < I < 21 degrees and \b\ < 1.5 degrees. The catalogue has 
a 99% completeness for sources with a peak flux above 6cr, which corresponds to a flux density of ~0.4Jybeairr'. The parameters 
extracted for sources with peak fluxes below the 6<x completeness threshold should be used with caution. Tests on simulated data 
find the uncertainty in the flux measurement to be ~12%, however, in more complex regions the flux values can be overestimated by 
a factor of 2 due to the additional background emission. Using a search radius of 30" we found that 40% of ATLASGAL compact 
sources are associated with an IRAS or MSX point source, but, ~50% are found to be associated with MSX 21 /jm fluxes above the 
local background level, which is probably a lower limit to the actual number of sources associated with star formation. 
Conclusions. Although infrared emission is found towards the majority of the clumps detected, this catalogue is still likely to include 
a significant number of clumps that are devoid of star formation activity and therefore excellent candidates for objects in the coldest, 
earliest stages of (high-mass) star formation. 

Key words. Stars: formation - Surveys - Submillimeter - Catalogues 



1. Introduction 

The processes and time scales involved in the formation of high- 
mass stars are still poorly understood, although such stars play 
an important role in the evolution of their parental clouds and the 
subse quent star formation inside them (iKrumholz & BonnellL 
l2007l) . Moreover, they interact with their environment injecting 
an enormous amount o f energy both radiative and mechanical 
(iHoare & Francoll2007l) . 

Great advance s have been made i n unde rstanding how low- 
mass stars form dMcKee & Ostrikeil l2007t) . but the situation 
for high-mass stars is still very unclear. Different scenarios 
have been proposed that provide a theoretical framework for 
the formation of high -mass stars such as competitive accretio n 
dBonnell et al.Ll2004l) and core accretion dMcKee & TanLl2003l) 



On the observational side, the short time scales jMcKee & TanL 



* The full catalogue and the calibrated emission maps 
are only available in electronic form at the CDS via 
anonymous ftp to cdsarc.u-strasbg.fr (130.79.125.5) or via 
http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/ 

** Current address: CSIRO Astronomy and Space Science, P.O. Box 
76, Epping NSW 1710, Australia, e-mail: yanett@das.uchile.cl 



120021) involved in the different stages of the formation process 
make it difficult to build large representative samples for the sta- 
tistical studies necessary to robustly test these theories. 

Several surveys have been carried out to study the environ- 
ments of high-mass star-forming regions. The drawback of many 
of these surveys is that they target samples that fulfil certain 
selection criteria (e.g. maser emission, or peculiar far-infrared 
colours), and, as suc h, they are biased towards given evolution- 
ary stages (see e.g. IWvrowskiL 120081) . For example, methanol 
masers have been found to be exclusively associated with the 
very early stages of high-mass star formation up to the more 
evolved ultracompact HII region phase an d were the subjec t 
of the recent methanol multibeam (MMB; Gre en et al. 120091) 
surve y. The Red MSX Source survey (RMS; lUrquhart et al.l 
120081) targets high-mass young stellar objects (MYSOs) that are 
known to be mid-infrared bright, while the Co-Ordinated Radio 
and Infrared Sourc e for High-Mass Star Formation (CORNISH; 
IHoare etaD 120121) targets the ultra-compact (UC) Hn region 
stage at which point the young embedded star begins to ionize its 
surrounding envelope and hence represents a later evolutionary 
stage. 
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All of these surveys provide useful statistics on the partic- 
ular stage targeted (an evolutionary snapshot) but do not pro- 
vide a complete overview of high-mass star formation. What is 
needed is a Galaxy wide survey using an unbiased tracer that 
links all of these complementary surveys together. The ideal 
tracer to choose is the submillimetre thermal emission from dust 
as it is optically thin, and so traces the whole column of material 
across the Galaxy, and is associated with all stages of high-mass 
star formation. However, until recently, limitations in mapping 
speeds have resulted in limited coverage of the Galactic plane. 
For example, SCUBAfl 9-year life time, a total area of o nly 29.3 
deg 2 has been surveyed at 850 /urn dDi FrancescoLl2008l) . 

The need for unbiased surveys to get a complete census of 
the regions where high-mass stars are being formed has led to 
several projects to map the Galactic plane at millimetre and 
submillimetre wavelengths. The AP E?fl Telescope Large Area 
Survey of the Galaxy (ATLASGAL: ISchuller etalll2009h is the 
first complete survey of the inner Galactic plane providing an 
unbiased view and complete census of high-mass star-forming 
cores and clumps. It is well co mplemented by the H erschel 
Infrared Galactic plane (Hi-GAL; iMolinari et al.L 1201 Oh survey 
in the far-infrared range at a similar high spatial resolution of 
~20" (e.g., 0.5 pc at 5kpc), and at slightly long er wavelengths 
with the Bolocam Galactic Plane Survey (BGPS: lA.guirre et al.L 
1201 lh . Sub-mm surveys, as opposed to other surveys at shorter 
wavelengths such as the Galacti c Legacy Infrared Mid -Plane 
Survey Extraordinaire (GLIMPSE IBeniamin etall 120031) or the 
MIPS Inner Galactic Plane Survey (MIPSGAL: ICarev et al.L 
I2009h . are sensitive to dust in a broad range of temperatures, 
including the coldest dust. Therefore, such surveys have the po- 
tential to trace the earliest stages of star formation. 

In this paper, we present a compact source catalogue 
extracted from a limited range of the ATLASGAL survey. 
ATLASGAL is an unbiased survey of the Galactic plane 
at 870 nm w i th LA BOCA (Large APEX Bolometer Camera; 
ISiringo et al.L 120091) . and is the largest and most sensitive 
ground-based sub-mm survey to date. The catalogue presented 
here is extracted from a subset of the data collected from 2007 
to 2010, which covered 153 deg 2 , about one third of the total sur- 
vey area (420 deg 2 ). A catalogue will be compiled from the re- 
maining fields and be published in a companion paper when the 
data reduction and analysis of these fields has been completed. 

The aim of this survey is to obtain a complete census of 
cold dusty clumps in the Galaxy, and to study their distribu- 
tion across the Galaxy. This survey reveals a complete sam- 
ple of different stages of high-mass star formation. It will be 
complem entary to the other G alactic plane surveys such as 
UKIDSS (lLawrence et al.l2007|) and VISTA Variables in the Via 
Lactea (VW lSaito et al.ll2012l) at the near-infr ared; GLIMPSE 
and M IP SGAL at the mi d- to far-infrared (IBeniamin et alJ 
120031 and ICarev et al.ll2009l respectively); Hi-GAL which spans 
the far-infrared and sub-millimetre wavelengths (70-500 /im; 
IMolinari et aT1l2010l) between MIPSGAL and ATLASGAL; and 
finally the 5 GHz CORNISH Survey. These surveys combined 
together will provide a complete wavelength coverage of the 
Galactic plane from the near-infrared to the radio, and will also 
form a unique database for follow-up observations. 



To demonstrate the potential impact of this unique database, 
the ATLASGAL da ta has been used in numerous studies; 
IWienen et al.l (1201 2|) determined the temperatures and kine- 
matic distances of ~1000 clumps in various evolutionary stages 
using ammonia follow-up observations; ~1500 ATLASGA L 
sources are part of the MALT90 project dFoster et al.L 1201 llP l 
which simultaneously maps these sources in 16 spectral lines 
in the 3 -mm waveleng t h rang e with the Mopra radio tele- 
scope; iDeharveng et al.l d2010l) studied the association of in- 
frared bubbles wi th ATLASGAL dust clu mps; and based on 
ATLASGAL data iTackenberg et al.l d2012l) identified a number 
of starless cores, using indicators of star formation in GLIMPSE 
and MIPSGAL images. 

In Sect. [2] we describe the observations, and in Sect.[3]we 
detail how the sources were extracted. Tests on the reliability of 
the extraction algorithm are described in Sect. [4] In Sect. [5] we 
present the compact source catalogue and discuss the distribu- 
tions of the various source parameters, while in Sect.[6]we inves- 
tigate the correlation of ATLASGAL sources with the IRAS and 
MSX source catalogues and with other types of objects. Finally, 
we summarise our results in Sect.|7] 



2. Observations 



1 S ubmillimeter Common User Bolometer Array, ^Holland et all 
Il999h 

2 This publication is based on data acquired with the Atacama 
Pathfinder Experiment (APEX) is a collaboration between the 
Max-Planck-Institut fur Radioastronomie, the European Southern 
Observatory, and the Onsala Space Observatory. 



The o bservations were carried out with LABOCA dSiringo et all 
12009b at the APEX ( Atacama Path f inder E xperiment) 12 m sub- 
millimetre antenna dGiisten et al.L 120061) . located in Llano de 
Chajnantor, Chile. The typical pointing r.m.s. error was mea- 
sured to be ~4". The bandpass of LABOCA is centred at 
345 GHz, with a bandwidth of 60 GHz. At this frequency, the 
APEX beam FWHM is \9'.'2. 

In the present paper, we will focus on the region 330° < 
I < 21°, -1.5° < b < 1.5° (153 deg 2 ). Future versions of 
the catalogue will cover the full 420 deg 2 observed as part of 
ATLASGAL. 

The observations were made using on-the-fly maps with a 
scanning velocity of 3' /s. They produced slices of 1° width in 
Galactic longitude times 2° length in Galactic latitude for the 
observations carried out during 2007 and slices of 1 in Galactic 
longitude times 2 to 3° in Galactic latitude during 2008 to 2010. 
To avoid artefacts inherent to the observing technique, each posi- 
tion in the sky is covered at least twice, in two or more maps ob- 
served with different scanning angles. The flux calibration was 
made using planets as primary calibrators (Mars, Neptune and 
Uranus), and secondary calibrators (bright Galactic sources with 
known fluxes from the commissioning of LABO CA). The errors 
in pea k flux are estimated to be lower than 15% (ISchuller et all 
120091) . 

The data reduction wa s made us i ng the BOlometer array data 
Analysis package (BoA; ISchullerL 120121) . The steps involved 
in the data reduction are explained in detail in ISchuller et al.l 
(120091) . and can be summarised as: flagging bad bolometers, re- 
moving the correlated noise, flagging noisy bolometers, despik- 
ing, low frequency filtering, and subtracting a first order base- 
line. The maps are then built by co-adding the signals from all 
bolometers, using a weighted mean with natural weighting. A 
weight map is also calculated from the sum of all the weights 
that contribute to each pixel in the final map. An r.m.s. map can 
then be computed as l/sqrt( weight). 

A problem inherent to the reduction of ground-based bolo- 
meter data is the loss of the extended emission, which is filtered 
out when the correlated noise is subtracted. As a result, uniform 



http://malt90.bu.edu/ 
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Fig. 1. In the left panel we present an r.m.s. noise map of a 3°x3° field centred on £ = 19.5° and b — showing how the noise varies 
as a function of Galactic longitude and latitude. In the right panel we present a histogram of the pixel values found within \b\ < 1° 
where the noise is approximately constant. The red line show the result of a Gaussian fit to the noise, from which we estimate the 
average noise for this field (cr rm s - 60mJy). 



emission at angular scales larger than 2'.5 is removed. This re- 
duction scheme was optimised to recover compact sources. The 
final maps are built with a pixel size of 6", corresponding to a 
~ 1/3 of the telescope beam. The size of each map is of 3 x 3 de- 
gree, having an overlap between two consecutive maps of ~4.5 
arc minutes. 

In the left panel of Fig.Q] we present an example of a fi- 
nal noise map, i. e. the r.m.s. noise level distribution over a re- 
duced map. As can be seen from this map the noise is reason- 
ably uniform within \b\ < 1°, however, due to poorer coverage 
at higher latitude (\b\ > 1°) the pixel to pixel r.m.s. is signifi- 
cantly larger, typically 150-250mJy beam -1 . The vast majority 
of the emission seen in the ATLASGAL maps is concentrated 
close to the Galactic mid-plane and so the sensitivity of the sur- 
vey can be considered fairly uniform. In the right panel of Fig.Q] 
we present a histogram of the flux values of pixels with \b\ < 1°. 
The r.m.s. noise is well modelled by a Gaussian fit (shown as a 
dashed red line), from which we estimate the r.m.s. noise to be 
~ 60mJy beam -1 . We estimate the r.m.s. noise in each field in 
the same way and present a plot showing how the noise varies as 
a function of Galactic longitude in Fig. [2] This plot reveals that 
the noise is fairly uniform for the regions presented here (~ 50- 
70 mJy beam -1 ), with the observed variations resulting from dif- 
ferent weather conditions and elevation of the various observa- 
tions. This is particular noticeable for the maps within +7 and 
+ 17 degrees where the higher noise in these maps is believed to 
be due a combination of high airmass and poorer weather condi- 
tion. 



90^" 
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Fig. 2. Plot showing how the noise varies as a function of 
Galactic longitude. The points show the standard deviation of 
the pixel to pixel variation averaged for each 3°x3° field, how- 
ever, we have only included pixels with \b\ < 1° where the noise 
is relatively uniform and where the vast majority of the 870 fim 
emission is concentrated. A dashed line connects the point to 
emphasis the field-by-field variations in the noise. 



3. Source extraction 

To create the compact source catal ogue we used the sourc e 
extraction algorithm SExtractor dBertin & Arnoutsl Il996h . 
Although SExtractor was originally created as a method of 
source extraction for visible and infrared images, it has been 
used w ith good results for the extraction of sub-mm sources as 
well bv lThompson e t al. ( 2006) and more recently by the Planck 
team to produce the higher frequency components of their com- 



pact s ource catalogue (2 17-857 GHz; IPlanck Collaboration et alj 
1201 11) . 

Most source extraction algorithms assume that the noise 
level is constant across the emission map supplied. However, 
in most surveys this is unlikely to be the case due to variations 
in coverage, weather, elevation and instrumental effects, which 
make the noise level vary across the field, and this is indeed the 
case for the ATLASGAL survey. In Fig. Q] we present the noise 
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Fig. 3. Example of an ATLASGAL emission map (left panel) and the corresponding filtered signal-to-noise ratio map (see text for 
details). The region outlined in red indicates the field presented in Fig. [4] 



map for the 3° x 3° field centred on I =19.5° and b = 0. The fact 
that most source extraction algorithms assume a constant noise 
can lead to real sources lying in areas of lower than average noise 
regions being missed, while conversely, in higher noise regions 
there is a danger that spurious noise spikes will be mistakenly 
interpreted as real sources. 

In order to reduce the number of spurious sources and avoid 
missing genuine sources we have converted each emission map 
into a signal-to-noise ratio map by dividing it by its correspond- 
ing noise map. These maps also contain a large fraction of low 
surface brightness diffuse and filamentary emission, particularly 
towards the Galactic centre and towards the tangent with the 
Norma spiral arm and other prominent star forming regions. This 
not only complicates the identification of compact sources but 
can also lead to the identification of a large number of spuri- 
ous sources. To limit the source confusion due to large scale ex- 
tended diffuse emission we have removed it using the program 
FINDBACK, which is part of the CUPIE0 clump identification and 
analysis package. This produces maps with a higher contrast be- 
tween compact sources and their local background. In Fig.[3]we 
present an emission map to show the affect of the varying noise 
level in the left panel, while in the right panel we show the same 
map after smoothing the variations by dividing through by the 
noise map and spatial filtering to remove the diffuse background. 
These two images nicely illustrate the advantages of using the 
signal-to-noise map for the source detection. 

To summarise, we first smoothed the noise variations by di- 
viding the maps by the noise map and created a signal-to-noise 
map, and as a second step we applied spatial filtering to remove 
part of the extended emission. This is compared to the original 
maps with varying noise level in Fig. [3] 

SExtractor can work on two images simultaneously, one 
for the source detection and one for measuring the physical 
parameters of each source. We have used the spatially filtered 
signal-to-noise map for source detection, however, we have used 
the original emission maps to extract the corresponding source 
parameters. In this way, the flux and size of the sources are not 



4 CUPID is part of the STARLINK 
(http://starlink.jach.hawaii.edu). 



software suite 



affected by the filtering process performed for the source de- 
tection. These steps should ensure that only genuine sources 
make it into the final catalogue. In Fig.|4]we present a subre- 
gion, taken from the emission map presented in the left panel of 
Fig-E] showing the positions and sizes of sources identified by 
SExtractor. 

The ATLASGAL maps were searched for emission above a 
threshold of 3<x r , ms ., where <x r m s is the local noise level deter- 
mined from a user defined box by SExtractor. The value of 
this threshold was determined empirically by varying its value 
for a number of test regions so as to optimise the number of gen- 
uine sources detected while avoiding the inclusion of the more 
dubious sources. In addition to the threshold we also reject any 
sources with fewer pixels than the beam integral (i.e., ~11 pix- 
els). 

To separate possible blended sources, SExtractor divided 
the emission into 40 sub-levels of flux. In each of these 40 lev- 
els, the ratio of the integrated intensities of adjacent pixels is 
computed. If the ratio between two pixels is greater than 0.001, 
then those pixels are considered part of a different source. If not, 
both pixels are considered to be part of the same source. This 
comparison is done for each pixel in each level. To illustrate this 
process we present a schematic diagram in Fig. [5] that shows a 
complex emission profile that is divided in 23 sub-levels, at the 
10th level the source is separated in two branches then at the 12th 
level the first branch is separated again in two sources A and B, 
at the level 17th there is a small separation at the branch B, but 
the ratio of the integrated intensity of this pseudo new branch is 
lower than the ratio defined to separate sources so it is declared 
as part of the B branch. The process continues until the last level 
at which point seven discrete sources have been identified. 

4. Artificial sources 

In order to check the reliability of the compact source catalogue, 
we test our extraction algorithm on maps with artificial sources 
with known properties. To test the influence of the background 
emission in the extraction of sources, we created maps with uni- 
form and varying background noise. These two cases allow us 
to compare the performance of the extraction algorithm in the 
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Fig. 4. Example of the sources extracted. The field shown in greyscale is the region outlined in red in the left panel of Fig. [3] 
The positions and sizes of sources identified by SExtractor are shown as red and blue ellipses. The red ellipses indicate to the 
most reliable detections, while the blue ellipses indicated the detections associated with higher flag values, given by SExtractor, 
which correspond to less reliable detections (see discussion in Sect. l5.1l for details). The ATLASGAL beam size is shown in green 
in the lower left corner of the map. Contours start at 3<x (i.e., 0.18 Jy beam 1 ) and increase with steps of 0.48, 1.61, 4.32, 8.96, 
15.89 JybearrT 1 . 



idealised case where the background is smooth and the more 
realistic case when the sources are intermixed with large scale 
diffuse emission. 

The varying background noise map was created from the 
residual maps produced by SExtractor. The field used in this 
test is the map located at 19.5°, which is fairly representative 
of the rest of the ATLASGAL maps. The residual map contains 
background emission and zeros at the position of the extracted 
sources. The test map was smoothened using a boxcar func- 
tion with a width of 10 pixels to remove the sharp ridges left 
by the source removal. These steps produced an emission map 
with varying noise properties similar to the emission seen in the 
ATLASGAL maps. 

We used the CUPID task MAKECLUMPS to generate a catalogue 
of artificial objects with a uniform distribution in t, b. The flux of 
the injected objects was uniformly distributed between 0.01 and 
0.75 Jy beam -1 , however, the source size was set to a FWHM of 
5 pixels (30") with a Gaussian shape. Taking these input param- 
eters and using the test field as a template, MAKECLUMPS allowed 
us to produce a field of artificial sources with some added noise 
that has the same dimensions and coordinates of the test field. To 
produce the varying background test field we did not add addi- 
tional noise with the algorithm but simply co-added the artificial 



source map with the varying background noise map described 
in the previous paragraph. For the uniform background test field 
we simply used the output produced by the MAKECLUMPS task 
but this time set the noise level to be similar to the original map 
(i.e., ~60mJy beam -1 ). 

We repeated this process to produce 20 artificial source 
maps, all with both a varying and uniform background noise 
level. This group of 20 maps produced a catalogue of ~5000 
artificial sources, which is comparable to the number of gen- 
uine sources detected in the region presented here. For the 
varying background noise maps we produced background sub- 
tracted signal-to-noise maps following the same procedure used 
for the real ATLASGAL maps (see Sect.[3]for further details). 
SExtractor was then run on these artificial source maps using 
the same parameters as those used for the construction of the 
real catalogue. To compare the dependence of completeness on 
the background, the same extraction algorithm was tested on the 
maps with uniform noise. However, for these there was no need 
to try and filter the large scale variations in the background and 
therefore SExtractor was used directly on the simulated maps. 

In Fig.|6]we plot the number of sources recovered as a func- 
tion of input peak flux for the uniform and varying background 
noise maps (indicated by the blue and red lines, respectively). 
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Fig. 5. Example of the deblending process. The dependence of 
a pixel to one source or another will depend on the relative in- 
tegrated intensities (grey area) of their adjacent pixel for each 
level (dotted lin es). This figure has been reproduced here with 
permission from lBertin & Arnoutsl(ll996l) . 
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Fig. 6. Detection ratio of recovered to injected sources as a func- 
tion of input peak fluxes. The blue solid and red dashed lines 
show the detection ratios for sources injected onto a uniform and 
varying background noise map, respectively. The dashed vertical 
lines indicate (from left to right) the lcr, 3<x, 5cr, and lcr noise 
levels. 



In the case of the uniform noise level the source extraction al- 
gorithm is more than 99% complete above ~ 0.3 mJy beam" 1 , 
which corresponds to approximately 5<x. The shape of the de- 
tection ratio and the 5<x completeness level found for the uni- 
form back ground noise maps is ve ry similar to that found for 
the BGPS dRosolowskv et al.ll2010l) . However, the completeness 
profile determined from the sources injected and recovered from 
the varying background noise, which is much more realistic, is 
very different. In the case of the uniform background noise maps 
no sources are detected below the ~ 3<x level, however, we start 
to detect sources at ~ lcr when adopting varying -background- 
noise maps. The detection of low peak flux sources is better in 
the varying background noise maps compared to those with uni- 
form noise. This is because the peak fluxes are boosted by the 
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Fig. 7. In the upper panel we present a comparison of the re- 
covered values for the peak fluxes values to the input values of 
the injected sources on the maps with varying background. The 
blue line represents equality. In the lower panel we present the 
recovered to injected flux ratio (black histogram) and the results 
of a Gaussian fit to values less than 1.2 (dashed red line). 



diffuse continuum emission they have been injected onto. As a 
consequence, the fluxes of weaker sources (<5cr) are somewhat 
less reliable in this case. 

We also note that for the varying background noise case 
the ~99% completeness is obtained for peak fluxes above 
~ 0.4 mJy beam -1 , which, at approximately 6<x, is slightly worse 
than that found for the uniform background noise case. This 
difference can be attributed to the presence of the diffuse low 
surface brightness emission which makes it more difficult to 
recover sources below <6<x. Clearly the presence of a varying 
background due to large scale diffuse emission has a significant 
impact on the detectability and reliability of fluxes for weaker 
sources (i.e., <6cr) and therefore we do not advise to use the 
sources below 6cr blindly. 

In the upper panel of Fig.|7]we compare peak flux of the re- 
covered and injected sources. It is clear that these flux values are 
very well correlated and the injected and recovered fluxes are in 
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excellent agreement with each other (the blue line plotted over 
the data shows indicates equality). We also note the presence 
of a number of sources that have significantly higher recovered 
fluxes than input fluxes, which is a result of flux boosting from 
the varying background mentioned earlier. In the lower panel of 
Fig. [7] we present a histogram of the ratio of the injected and 
extracted flux values for each source. The distribution profile 
below a ratio value of 1 .2 can be approximated by a Gaussian 
fit to the data (shown as a dashed red line on the plot) with a 
peak at ~ 1 and has a <x of 0.12, which would suggests that the 
flux values extracted for the sources in this catalogue are gener- 
ally well determined with an uncertainty of 12%, however, for a 
small number of sources in more complex regions the flux can 
be overestimated by up to a factor of 2. 

5. Compact source catalogue 

5.1. Catalogue description 

SExtractor was run on each of the seventeen 3° x 3° fields 
that make up the inner most portion of the ATLASGAL survey 
(i.e., 330° < 6 < 21°). The results obtained for each field were 
subsequently combined in order to compile a complete catalogue 
of sources, taking care to eliminate duplicated sources located in 
the overlap regions between adjoining fields. In total we have 
identified 6639 sources in the 153 deg 2 . 

In TableQ] we present a sample of the catalogue (for the re- 
gion presented in Fig.[4j> as an example of its form and content. 
The full and up-to-date catalogue and ATLASGAL emission 
maps will be hosted by the ATLASGAL project pagfl the ESC0 
and CDS0 database. There will be incremental improvements in 
image processing, calibration and source extraction. In the ta- 
ble we present the source name, the source peak and barycentric 
position, the size of the semi-major and minor axis and their 
position angle, the deconvolved radius, the peak and integrated 
fluxes and their associated errors and any warning flags gener- 
ated by the algorithm (see full bit per bit description below). 

The source names are based on the Galactic coordinates of 
the peak flux position, which are given in Cols. 2 and 3. The 
barycentric coordinates are determined from the first order mo- 
ments of the longitude and latitude profiles. The barycentric po- 
sition generally defines the centre position of a source, however, 
if the emission profile of a source is strongly skewed or is associ- 
ated with substructures, the peak position can be a more accurate 
measure of a source's position. The source sizes describe the de- 
tected source as an elliptical shape, although no fitting is actually 
performed. The semi-major and semi-minor axis lengths repre- 
sent the standard deviation of the pixel co-ordinate values about 
the centroid position, weighted by the pixel values. The position 
angle is given as anti-clockwise from Galactic north. 

Since the source sizes are determined only using the pixels 
above the detection threshold it is likely these underestimate the 
true source sizes, and in some cases r esults in source sizes that 
are smaller than the beam. Following iRosolowskv et all d2010l) 
we estimate the angular radius (8r) from the geometric mean 
of the deconvolved major and minor axes and multiplied by a 
factor n that relates the r.m.s. size of the em ission distribution 
of the source to its angular radius (Eqn.6 of IRosolowskv et aTl 
12011 : 

0R = 4« j -oL)(c4 in -£^ m )l 1/4 , (1) 

5 www.mpifr-borvn.mpg .de/div/atlasgal/ 

6 Available through http : //archive . eso . org/wdb/wdb/adp/phs 

7 |http : //cdsweb . u- strasbg . fr/cgi-bin/qcat? J/A+A/| 



where cr fe „ is the r.m.s. size of the beam (i.e., <x/, m = 
S'fwhm/ V8 In 2 a 8")- The value of r/ depends on the size of 
the s ource with respect to the beam and the emission distribu- 
tion. IRosolowskv et al.l d2010t) adopted a value of 2.4, which is 
the median of the values derived from a range of models, how- 
ever, they note that the value of r\ can vary by a factor of 2 in 
their simulations. For consistency and to facilitate comparisons 
between the BGPS and ATLASGAL catalogues we also adopt 
this value, which is approximately equivalent to the sources ef- 
fective radi us (i.e., R e fi = ^/jA/n ), where A is the surface area of 
the source: iDunham et al.ll201 lb . 

The peak flux is directly obtained from SExtractor, how- 
ever, since the algorithm treats each pixel in the given emission 
map in a statistically independent manner and the ATLASGAL 
maps are gridded with ~3 pixels per beam it is necessary to ac- 
count for the beam area to determine the final measured fluxes 
and their associated uncertainties. This was simply done by di- 
viding the derived flux values by the beam area (i.e., 1.133 x 
FWHM 2 11.6 pixels). The errors given for the peak and inte- 
grated fluxes include the absol ute calibration uncerta inty of 15% 
mentioned in Sect. 2 (see also ISchuller et all 120091 for details) 
combined in quadrature with the intrinsic measurement error. In 
the case of the integrated flux this is provided by SExtractor, 
however, the error in the peak flux is derived from the local r.m.s. 
noise in the image. 

In Col. 15 we present the quality flag given by SExtractor 
which gathers all the infor mation about possible pr oblems or 
artefacts affecting the source dBertin & Arnoutslll996l) . For each 
source, this flag is either zero (no particular problem), or is equal 
to the sum of one or more number(s) with the following mean- 
ings: (1) the object has neighbours, bright and close enough to 
significantly bias the photometry, or bad pixels (more than 10% 
of the integrated area affected); (2) the object was originally 
blended with another one; (4) at least one pixel of the object is 
saturated (or very close to); (8) the object is truncated (too close 
to an image boundary); (16) object's aperture data are incom- 
plete or corrupted; (32) object's isophotal data are incomplete 
or corrupted; (64) a memory overflow occurred during deblend- 
ing; (128) a memory overflow occurred during extraction. For 
example, a flag value of 10 means that the object was originally 
blended with another source and that it is truncated because is 
located too close to the edge of the map. In the final column of 
the table we present the signal to noise for each source. 

In Fig. [4] we overplot the positions of sources detected by 
SExtractor towards a small region to demonstrate its perfor- 
mance. To check the reliability of sources with high flag values 
we plot sources with flag values less than 4 (which are the most 
reliable) in red while those with higher flag values are shown in 
blue. In this example region it is clear that the sources with the 
higher flag values are associated with genuine emission, how- 
ever, the emission does appear to be weaker and more diffuse 
than found to be associated with the lower flag value detections. 
A similar situation is seen in other regions examined and so 
we conclude that these high flag detections do identify genuine 
sources of emission, but with the caveat that their associated pa- 
rameters are somewhat less reliable. 

5.2. Catalogue properties 

In this section we will look at the Galactic distribution of sources 
detected and the overall distribution of their derived parameters. 
eAsmaid>e£rkiriB n the reliability of the ATLASGAL catalogue we 
compare our results to those obtained from the Bolocam Galactic 
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Table 1. The ATLASGAL compact source catalogue. The columns are as follows: (1) name derived from Galactic coordinates of 
the maximum intensity in the source; (2)-(3) Galactic coordinates of maximum intensity in the catalogue source; (4)-(5) Galactic 
coordinates of emission centroid; (6)-(8) semi-major and semi-minor size and source position angle measured anti-clockwise from 
Galactic north; (9) effective radius of source; (10)-(13) peak and integrated flux densities and their associated uncertainties; (14) 
SExtractor detection flag (see text for details); (15) signal to noise ratio (SNR) — values for sources with peak flux below 6<x 
detection should not be used blindly. 



Name 




^max 


^max 


t 


b 


v maj 




PA 


Sr 


" peak 


AS peak 


S int 


AS int 


Flag 


SNR 






0) 


(°) 


C) 


(') 


(") 


(") 


(°) 


(") 


(Jybeam ') 


(Jy) 








(1) 




(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


AGAL019.486- 


-00.199 


19.486 


-0.199 


19.487 


-0.198 


27 


17 


73 


48 


0.80 


0.13 


9.35 


1.67 


2 


12.9 


AGAL019.508- 


-00.447 


19.508 


-0.447 


19.509 


-0.447 


22 


13 


161 


34 


0.48 


0.10 


4.06 


0.86 





7.6 


AGAL019.526- 


-00.086 


19.526 


-0.086 


19.526 


-0.089 


1 1 


7 


14 




0.35 


0.08 


1.59 


0.45 





5.6 


AGAL019.539- 


-00.456 


19.539 


-0.456 


19.538 


-0.458 


26 


14 


87 


40 


0.76 


0.13 


7.72 


1.42 


2 


11.2 


AGAL019.572- 


-00.101 


19.572 


-0.101 


19.574 


-0.100 


18 


14 


4 


33 


0.71 


0.13 


4.13 


0.87 


2 


8.7 


AGAL019.581- 


-00.069 


19.581 


-0.069 


19.582 


-0.068 


13 


8 


176 


1 1 


0.48 


0.10 


3.90 


0.83 


3 


6.5 


AGAL019.582- 


-00.091 


19.582 


-0.091 


19.582 


-0.091 


12 


6 


142 




0.59 


0.13 


1.81 


0.49 


3 


6.5 


AGAL019.589- 


-00.079 


19.589 


-0.079 


19.589 


-0.079 


18 


12 


150 


28 


0.73 


0.14 


9.27 


1.66 


3 


8.3 


AGAL019.601- 


-00.092 


19.601 


-0.092 


19.601 


-0.094 


8 


6 


24 




0.40 


0.11 


1.49 


0.43 





4.2 


AG AL01 9.606- 


-00.072 


19.606 


-0.072 


19.607 


-0.074 


17 


12 


22 


28 


0.46 


0.11 


6.35 


1.22 


3 


5.6 


AGAL019.609- 


-00.137 


19.609 


-0.137 


19.609 


-0.138 


15 


1 1 


169 


23 


2.79 


0.43 


7.59 


1.40 


3 


32.3 


AGAL019.609- 


-00.234 


19.609 


-0.234 


19.609 


-0.234 


22 


19 


151 


45 


19.20 


2.88 


44.85 


7.01 


3 


252.3 


AGAL019.611- 


-00.121 


19.611 


-0.121 


19.611 


-0.120 


13 


10 


147 


18 


1.17 


0.20 


2.86 


0.66 


3 


13.8 


AGAL019.612- 


-00.259 


19.612 


-0.259 


19.612 


-0.258 


28 


16 


153 


46 


1.88 


0.29 


14.88 


2.51 


3 


25.7 


AGAL019.628- 


-00.264 


19.628 


-0.264 


19.626 


-0.264 


14 


7 


60 




0.46 


0.10 


1.62 


0.45 


3 


6.8 


AGAL019.631- 


-00.162 


19.631 


-0.162 


19.630 


-0.162 


12 


9 


175 


13 


0.70 


0.14 


5.80 


1.13 


3 


7.9 


AGAL019.641- 


-00.172 


19.641 


-0.172 


19.641 


-0.172 


19 


13 


32 


32 


0.53 


0.12 


11.06 


1.93 


2 


5.7 


AGAL019.686- 


-00.127 


19.686 


-0.127 


19.684 


-0.129 


25 


14 


103 


40 


0.70 


0.14 


7.52 


1.39 


2 


8.2 


AG AL01 9.694- 


-00.174 


19.694 


-0.174 


19.694 


-0.175 


10 


8 


62 




0.44 


0.10 


1.02 


0.34 





6.0 


AGAL019.699- 


-00.266 


19.699 


-0.266 


19.699 


-0.262 


26 


21 


138 


52 


1.97 


0.30 


12.69 


2.18 


3 


27.7 


AGAL019.706- 


-00.239 


19.706 


-0.239 


19.706 


-0.239 


27 


14 


136 


41 


0.88 


0.15 


9.56 


1.71 


3 


12.1 


AGAL019.716- 


-00.171 


19.716 


-0.171 


19.717 


-0.172 


12 


8 


58 


10 


0.44 


0.10 


2.62 


0.63 





5.7 


AGAL019.726- 


-00.114 


19.726 


-0.114 


19.726 


-0.114 


14 


1 1 


10 


23 


1.09 


0.18 


7.61 


1.41 





12.5 


AGAL019.746- 


-00.171 


19.746 


-0.171 


19.745 


-0.172 


18 


14 


88 


33 


0.40 


0.09 


3.93 


0.84 


2 


6.0 


AG AL01 9.754- 


-00.129 


19.754 


-0.129 


19.755 


-0.128 


22 


1 1 


116 


30 


1.87 


0.29 


11.23 


1.96 


2 


20.9 


AG AL01 9.771 - 


-00.114 


19.771 


-0.114 


19.770 


-0.113 


22 


7 


45 




0.54 


0.10 


2.78 


0.65 


2 


8.3 


AGAL019.811- 


-00.342 


19.811 


-0.342 


19.812 


-0.340 


16 


9 


95 


15 


0.37 


0.08 


3.27 


0.73 


19 


6.0 


AGAL019.817- 


-00.516 


19.817 


-0.516 


19.816 


-0.515 


8 


6 


68 




0.33 


0.08 


1.20 


0.38 





5.3 


AGAL019.829- 


-00.329 


19.829 


-0.329 


19.830 


-0.331 


22 


18 


20 


43 


1.15 


0.19 


9.90 


1.76 


2 


17.2 


AGAL019.874- 


-00.092 


19.874 


-0.092 


19.874 


-0.092 


10 


6 


13 




0.32 


0.08 


1.08 


0.36 





5.3 


AGAL019.882- 


-00.534 


19.882 


-0.534 


19.885 


-0.534 


30 


17 


178 


50 


7.43 


1.12 


27.37 


4.39 


2 


90.7 


AGAL019.899- 


-00.441 


19.899 


-0.441 


19.896 


-0.439 


11 


9 


56 


10 


0.35 


0.08 


2.51 


0.61 


18 


5.7 


AGAL019.902- 


-00.582 


19.902 


-0.582 


19.902 


-0.582 


13 


7 


35 




0.38 


0.09 


3.24 


0.73 


16 


5.9 


AGAL019.906- 


-00.241 


19.906 


-0.241 


19.906 


-0.238 


13 


10 


44 


20 


0.45 


0.10 


3.33 


0.74 





5.8 


AGAL019.911- 


-00.206 


19.911 


-0.206 


19.911 


-0.205 


11 


10 


45 


15 


0.33 


0.08 


2.30 


0.57 





4.9 


AGAL019.922- 


-00.259 


19.922 


-0.259 


19.926 


-0.258 


25 


15 


103 


42 


2.51 


0.38 


13.96 


2.37 





36.8 


AGAL019.967- 


-00.106 


19.967 


-0.106 


19.970 


-0.105 


11 


9 


116 


14 


0.45 


0.09 


2.29 


0.57 





7.2 


AG ALO 19.977- 


-00.214 


19.977 


-0.214 


19.976 


-0.216 


16 


10 


158 


21 


0.87 


0.15 


3.12 


0.70 





13.9 


AGAL019.991- 


-00.169 


19.991 


-0.169 


19.991 


-0.169 


10 


7 


86 




0.55 


0.11 


1.52 


0.43 


3 


7.9 


AGAL020.004- 


-00.167 


20.004 


-0.167 


20.002 


-0.168 


13 


10 


115 


19 


0.40 


0.09 


4.63 


0.95 


18 


5.7 


AGAL020.019- 


-00.151 


20.019 


-0.151 


20.018 


-0.152 


9 


6 


119 




0.37 


0.09 


1.96 


0.51 


16 


5.0 


AGAL020.022- 


-00.169 


20.022 


-0.169 


20.024 


-0.170 


18 


7 


80 




0.38 


0.09 


5.80 


1.13 


19 


5.5 


AGAL020.081- 


-00.136 


20.081 


-0.136 


20.080 


-0.136 


18 


16 


106 


36 


7.68 


1.15 


18.43 


3.04 





119.8 



Notes: Only a small portion of the data is provided here, the full table is only available in electronic form at the CDS via anonymous ftp to 
cdsarc.u-strasbg.fr (130.79.125.5) or via http://cdsweb.u-strasbg.fr/cgi-bin/qcat7J/A-l-A/ 



Plane Sur vey (BGPS:lAguirre et al.l201 ll) that were recently pre- 
sented bv lRosolowskv et alj d2010l) . 

The BGPS covers 150 deg 2 of the Galactic plane, includ- 
ing the majority of the first quadrant with a latitude range of 
\b\ < 0.5°. Although the overlap between the BGPS and the 
ATLASGAL region presented in this paper is relatively mod- 
est (~30deg 2 ) when compared to the total area covered by the 
two surveys, the overlapping region does contain some of the 
highest density regions in the plane, and thus, provides a large 
number of sources in common to both surveys. Another thing to 



bear in mind when comparing the two catalogues is that the two 
surveys were preformed at different wavelengths and hence they 
have different spatial resolution and sensitivity and so we do not 
expect a one to one correlation between individual sources or 
their parameters. However, we do expect the overall distribution 
of source parameters to be correlated. To facilitate the compar- 
isons between the two catalogues we have followed the structure 
and present the same param eters in Table 1 as are gi ven in the 
BGPS catalogue (cf. Table k lRosolowskv et al.ll2010l) . 
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Fig. 8. Results of the cross-matching preformed between the ATLASGAL, BGPS and MMB catalogues. In the upper left panel we 
plot the peak fluxes for BGPS and ATLASGAL sources found within a matching radius of 40". In the upper right panel we compare 
the angular offsets in Galactic longitude and latitude between the ATLASGAL and BGPS catalogues. In the two lower panels we 
present plots showing the angular offsets in t and b between matched MMB-ATLASGAL (left) and MMB-BGPS (right) sources. 
The dashed lines overplotted in the in the upper right and lower panels show the zero offset positions in I and b, while the blue 
crosses indicate the mean offset of the matched samples. 



5.2.1 . Astrometry 

In order to check the astrometry of the ATLASGAL catalogue 
we compared the peak positions determined in our catalogue 
with those determined by the BGPS. Using an angular search 
radius of 40" we identify ~3, 000 likely associations. The match- 
ing radius was chosen to be a little larger than the resolution of 
the BGPS to obtain a large sample of matched sources to com- 
pare the positions, while it is still small enough to avoid multiple 
matches between sources. In the upper left panel of Fig. [8] we 
plot the peak 870 fim and 1 . 1 mm flux densities of the associ- 
ated sources. As this plot shows, there is a strong correlation 



between the fluxes of the matched sources which suggests that 
the matched sources are really associated with each other. 

In the upper right panel of Fig. [8] we plot the angular po- 
sition offsets in Galactic Longitude and Latitude between the 
matched ATLASGAL and BGPS sources. Although we should 
expect a spread in the angular offsets for individual matched 
pairs of sources, given the large number of matches we also ex- 
pect the mean of the offset distribution to be centred around zero. 
However, visual inspection of this plot reveals a systematic off- 
set between the two catalogues of ~5" in latitude (see Table|2]for 
details). This offset is comparable to the pixel sizes for the two 
surveys (6" and 7" for ATLASGAL and BGPS, respectively). 
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Table 2. Angular offsets between peak flux positions deter- 
mined from the ATLASGAL, BGPS and MMB catalogues. The 
errors given are the standard errors. 



Catalogues 


M 




Ab 




Matches 


ATLASGAL-BGPS 


-4.74 ± 


0.26 


+ 1.23 ± 


0.23 


2803 


MMB-ATLASGAL 


-0.20 ± 


0.51 


+0.79 + 


0.43 


459 


MMB-BGPS 


-4.30 ± 


0.86 


+0.92 ± 


0.70 


209 


SCUBA-BGPS" 


-1.8 ± 


1.2 


+0.4 ± 


0.8 





" These results have been taken from A guirre et all d2011h . 



Given that this offset is significantly smaller than their respec- 
tive beam size, and that the BGPS and ATLASGAL sources are 
not point sources, means that this systematic offset cannot be 
seen from a comparison of contoured emission maps of individ- 
ual sources. However, the offset is approximately 20 times the 
standard error and is therefore significant. 

In an effort to track down the origin of this offset we re- 
peated the matching procedure described at the beginning of 
this subse ction with the Me thanol Multibeam (MMB) survey 
catalogue ( Green et al 1 120091) . The MMB catalogue consists of 
~700 methanol masers observed at high resolution with the 
Australia Telescope Compact Array which has provided posi- 
tions accurate to better than an arcsecond r.m.s. (ICaswell et al.l 
l2010h . Methanol masers are associated with the early stages of 
high-mass star formation, which takes place within deeply em- 
bedded dust clumps. We should therefore expect a strong posi- 
tional correlation between the methanol masers and dust clumps 
identified by ATLASGAL. We have matched ~500 masers with 
ATLASGAL sources and present a plot of their positional off- 
sets in the lower left panel of Fig. [8] It is clear from this plot and 
from the mean offsets and standard errors presented in Table|2] 
that there is no significant offset between the two catalogues. 

We also computed the positional offsets between the BGPS 
and MMB catalogues (see lower right panel of Fig. [8). We find a 
very similar systematic offset (A( *5") as found from the com- 
parison of the ATLASGAL and BGPS catalogues shown in the 
upper right panel of Fig. [8] Although the number of matches 
found between the MMB and BGPS samples is smaller (-200) 
the mean offset in longitude is five times the standard error and 
is therefore significant. 

lAguirre etal] d2011l) checked the astrometry of the BGPS 
catalo gue against the SCUBA Legacy Catalogue dDi Francesco! 
12008b hereafter SLC) and found the offsets between them to be 
M = -1.8"±1.2" and Ab = 0.4"+0.8" J!|The mean offset va lues 
and their associated error reported by lAguirre et alj d201 lb are 
consistent with the zero offset position in longitude and latitude. 
However, because their standard errors are quite large their re- 
sults are consistent with the values determined here for the off- 
set between the ATLASGAL and BGPS sources (i.e., within 3 
standard errors). The excellent angular correlation between the 
MMB and ATLASGAL sources, and the fact that the same sys- 
tematic offset is observed when comparing the ATLASGAL and 
MMB catalogues with the BGPS, would suggest that the ob- 
served systematic offset originates from the BGPS catalogue. 



8 In the lAguirre et al .1 d2011h paper they give the offsets determined 
from the difference between the BGPS and SCUBA fields, however, to 
be consistent with the way the offsets have been determined here these 
values have been reversed. 



5.2.2. Galactic distribution 

In Fig. [9] we present histograms showing the Galactic longitude 
and latitude distribution of the ATLASGAL catalogue (left and 
right panels, respectively). Inspection of the longitude distribu- 
tion reveals a number of peaks in the source counts; these are lo- 
cated at approximately -27°, -22°, -8°, 0° with a weaker peak 
at 14°. The Galactic structure trace d by ATLASGAL so urces 
was the subject of a recent paper by iBeuther et al.l d2012l) . The 
authors used the CLUMPFIND source extraction algorithm to pro- 
duce a source catalogue and they analysed the distribution in 
Galactic longitude. They found that the positions of the peaks 
in source counts coincide with the positions of the tangent 
of spiral arms and individual star forming regions. Identifying 
the peaks in the source counts with the tangent of the Norma 
arm (-28°), the farside of the Near 3-kpc arm (-22°), the 
NGC6334/NGC6357 star forming regions (-8°), the Galactic 
centre region and finally M16 and M17 regions (+15°), respec- 
tively. 

Comparing the longitude distribution of ATLASGAL 
sources with the BGPS distribution we see that they are sim- 
ilar, both showing the same peaks, however, we note that the 
ratio of the Galactic centre peak to the other peaks is quite dif- 
ferent. This is likely due to the different spatial resolution of the 
two surveys. The complex emission seen towards the Galactic 
centre is more easily deblended in the ATLASGAL data due to 
the higher spatial resolution which results in a higher number 
of sources detected in this region. Conversely, the lower resolu- 
tion BGPS survey is more sensitive to lower surface brightness 
structures resulting in a higher proportion of sources being found 
away from the Galactic centre. 

The latitude distribution shown in the right panel of Fig.|9] 
peaks significantly below 0° as also noted by ISchuller et al.l 
d2009h and more recently bv lBeuther et al.l (120121) . This offset is 
consistent with the location o f the Sun above the Galactic plane 
(iHumphrevs & Lars en. 1995). Comparing the latitude distribu- 
tion to that of the BGPS we find them again to be similar. Both 
ATLASGAL and BGPS distributions have the same FWHM. We 
note that the tails of the distributions appear to be significantly 
different, however, the BGPS latitude range is \b\ < 0.5° with 
only a few excursions to larger value of \b\ and so the distribu- 
tions are not really comparable for latitudes \b\ > 0.5°. 

5.2.3. Flux distribution 

In Fig.[10]we show the peak and integrated flux distributions for 
the ATLASGAL sources (grey filled histogram) and the BGPS 
sources located in the overlap region (blue histogram). These 
distributions show the di fferential flux density sp ectrum (i.e., 
dN/dSv - AN/ASv oc S" a : lRosolowskv et al.ll2010l) . Comparing 
the distributions of the peak and integrated fluxes for the two cat- 
alogues they seem to be very similar, however, with the higher 
frequency measurements being shifted to higher values. 

The flux distributions above the turnover can be approxi- 
mated by a power-law that extends over three orders of mag- 
nitude. Fitting these parts of the distributions with a linear least- 
squares fit we find that the peak and integrated distributions can 
be represented by the values of a = 2.4 + 0.04 and 2.3 + 0.06, 
respectively. The results of these fits are indicated by the red 
lines overplotted on the left and right panels of Fig. [10] The val- 
ues obtained from the same fit to the BGPS p eak and integrated 
flux d istributions are 2.4 and 1 .9, respectively dRosolowskv et al.l 
l2010h . We note that the fit to the integrated fluxes is significantly 
steeper for the ATLASGAL catalogue than determined for the 
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Fig. 9. Histograms of the Galactic longitude and latitude distribution for ATLASGAL sources extracted by SExtractor (grey filled 
histogram) and the BGPS (blue histogram). In both plots the peak of the BGPS distribution has been normalised to the peak of the 
ATLASGAL distribution. The bin sizes used for the longitude (left panel) and latitude (right panel) distributions is 1° and 0.1°, 
respectively. 
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Fig. 10. Flux density number distribution for ATLASGAL sources extracted by SExtractor (grey filled histogram) and the BGPS 
(blue histogram). In both plots the peak of the BGPS distribution has been normalised to the peak of the ATLASGAL distribution. 
In the left and right panels we present histograms of the peak and integrated flux densities measured for each source respectively. 
The red line on each plot shows the result of a linear least-squares fit to flux densities larger than the peak of each distribution. 



BGPS catalogue. This probably results from the different reso- 
lutions between with two surveys as the lower resolution BGPS 
survey is likely to blend smaller clumps, that are resolved in 
ATLASGAL, together into larger clumps with higher integrated 
fluxes. This leads to a higher proportion of the BGPS sources 
being found at higher fluxes and the shallower slope than found 
for the ATLASGAL catalogue. 

The difference in the slope between the peak and integrated 
flux densities results from the fact that point sources dominate 
the peak distribution towards the completeness limit, whereas 
brighter sources tend to be larger and thus have larger total 
flux densities than for th e smaller sources. As pointed out by 
iRosolowskv et alJ (120101) . this makes the integrated flux density 



distribution top-heavy and effectively lowers the magnitude of 
the exponent. 

Another interesting feature of the peak distribution is that 
the main body of the distribution drops off completely at 
~70Jy beam -1 , however, there is another outlying peak located 
at ~150 Jy beam -1 . The bright sources that contribute to this flux 
density bin are the brightest 870 /urn sources in the catalogue and 
form part of a dust ridge associated with the Sgr B2 high-mass 
star forming region. 

5.2.4. Angular size distribution 

In Fig.QT]we present plots showing the size distributions of the 
catalogued sources. In the upper panel of this figure we show 
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histograms of the semi-major (cr ma j ) and semi-minor (cr min ) dis- 
tributions and the ATLASGAL beam size (<x bm 8", red dotted 
line). We note that there is a significant number of sources with 
sizes smaller than the beam, particularly with regard to the semi- 
minor axis. This is due to the fact that the sizes are derived from 
moment measurements of the emission distribution above the 
threshold value (~3cr rm s ) and therefore do not take account of 
all of the emission associated with eac h source down to the zero- 
intensity level. As previously noted bv lRosolowskv et ail d2010l) 
this has a disproportionate effect on weaker sources and often re- 
sults in their sizes being underestimated. We find the mean and 
median values of the peak fluxes for sources in our catalogue that 
are smaller than the beam are 0.52 and 0.43 Jy beam -1 , respec- 
tively, and therefore this sample is dominated by sources around 
the sensitivity limit with only a few percent having peak fluxes 
above 1 Jybeam -1 . 

In the middle panel of Fig.QT|we present a histogram show- 
ing the elongation distribution, that is the ratio of semi-major 
to semi-minor axes. The distribution has a peak and median 
value of 1.3 and 1.6, respectively, which indicates that most 
sources in the catalogue are significantly elongated. The elon- 
gation distribution closely matches that found for the BGPS 
sources. The only noticeable difference is that the ATLASGAL 
catalogue tends to have a slightly higher degree of elongation, 
possibly a result of the higher resolution of ATLASGAL survey. 
However, both catalogues appear to contain a significant number 
of very elongated sources reflecting the fundamentally filamen- 
tary structure of molecular clouds. 

In the lower panel of FigQj] we present the effective radii 
of the ATLASGAL and BGPS detections. The two distributions 
again appear to be broadly similar, however, the whole BGPS 
source distribution appears to have been shifted 10-15" to the 
right of the ATLASGAL distribution, which is clearly a result of 
the different survey resolutions. 



6. Comparison with other Galactic plane surveys 

6.1. Infrared catalogues 

We performed a cross correlation of o ur catalogue with the far - 
infrared IRAS Point Source Catalogue dJoint Iras Sciencelll994l) . 
and with the mid-infrare d Midcourse Spa ce experiment (MSX) 
Point Source Catalogue dEgan et al.ll2.Q03h . We used a search ra- 
dius of 30" around the peak flux position for all these survey 
catalogues. The resolution of the MSX survey and the APEX 
telescope used for ATLASGAL are very similar (18" and 19.2", 
respectively) making the correlation between both data sets rela- 
tively straightforward. This represents an initial attempt to asso- 
ciate ATLASGAL sources with their IR counterpart, and further 
associations with other catalogues (e.g. those extracted from the 
Spitzer GLIMPSE and MIPSGAL surveys) will be addressed in 
a future paper. 

Infrared emission detected by the IRAS and MSX surveys 
generally implies that star formation is already underway within 
a clump, therefore, these surveys trace objects in a relatively 
evolved stage of star formation. The ATLASGAL catalogue pre- 
sented here contains 6639 sources, of which we found 755 as- 
sociated with IRAS sources and 2314 associated with MSX point 
sources. The number of sources with IRAS and MSX association 
is 645. 

We found that 4216 objects, corresponding to 62% of the 
ATLASGAL sources, are not associated with any mid- or far- 
infrared point sources found in the MSX or IRAS catalogues. 
These objects are excellent candidates for cold dense clumps in 
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Fig. 11. Histograms of the angular size distribution of detected 
ATLASGAL sources. In the upper panel we show the distribu- 
tion of the semi-major (filled grey histogram) and semi-minor 
(blue open histogram) size distributions of detected ATLASGAL 
sources. The red dashed line shows the ATLASGAL r.m.s. beam 
size . The bin size used for both parameters is 4". In the middle 
and lower panels we present histograms of the elongation and ef- 
fective radii of ATLASGAL sources (grey filled histogram) and 
the BGPS (blue histogram). In both cases the peak of the BGPS 
distribution has been normalised to the peak of the ATLASGAL 
distribution. The bin sizes used for the elongation and radii dis- 
tributions are 0. 1 and 4" respectively. 
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Fig. 12. Distribution of separation between methanol maser 
sources from the MMB survey and ATLASGAL sources. This 
plot has been truncated at an angular offset of 50", however, 
only 7 ATLASGAL-MMB associations have larger offsets than 
this. The bin size is 2". 

an early stage of evolution, prior to the birth of high-mass pro- 
tostars. Further correlation of these sources with other infrared 
surveys such as GLIMPSE or MIPSGAL, can give further in- 
formation on their evolutionary stage, based on the presence or 
absence of star formation indicators such as green fuzzies (en- 
han cement at 4.5 //m, often re ferred to as extended green objects, 
e.g. JCvganowski et al.ll2008h or 24 /urn point sources. Moreover, 
follow up studies of these cold dense clumps in molecular lines 
and their detailed study with next generation telescopes such as 
the ALMA and the Jansky VLA will help improve our under- 
standing the initial stages of high-mass star formation. 

6.2. Methanol masers 

Methanol masers are well-known probes of early phases of high- 
mass s tar formation, in particular the strong 6.7 GHz Class-II 
maser (lMentenlll99ll) . Recently a Galactic plane survey of this 
maser transition usin g a 7-beam receive r on the Parkes telescope 
has been completed dGreen et al.L [2009). These maser detections 
were followed up at high angular-resolution with the Australian 
Telescope Compact Array to obtain sub-arcsecond positional 
accuracy. To date the MMB has reported the positions of 707 
me thanol masers located between 186° < I < 20° and \b\ < 
2° (ICaswell et al|l2010t IGreen et al.ll^oTok ICaswell et al.ll201 it 
iGreen et al.ll2012l) . of these 493 are located within ATLASGAL 
fields presented in this paper (i.e., 330°< I < 21°). Cross- 
matching these MMB sources with the ATLASGAL catalogue 
generated from the SExtr actor algorithm we find 471 posi- 
tional associations, which corresponds to ~96% of the available 
MMB sources but only -7% of the ATLASGAL catalogue. 

In Fig. [12] we show the angular separation between the 
ATLASGAL-MMB associations. The sharp peak seen towards 
small angular offsets reveals a strong correlation between the 
position of the methanol maser and the peak column density of 
each clump as traced by the submm continuum emission. We 
find that -87% of all ATLASGAL-MMB associations are found 
with angular separations < 12", after which the distribution of 
the offsets flattens to an almost constant background level. We 
find that 61 ATLASGAL-MMB associations have angular off- 



sets larger than 12"; this could indicate the presence of clumpy 
substructure that is smaller than the size of the APEX beam or 
a chance alignment with an MMB source associated with an un- 
detected dust clump. 

Inspection of the ATLASGAL maps at the locations of 
the 22 MMB sources not found to be associated with an 
ATLASGAL source reveals that nearly all of these masers are 
associated with either weak diffuse emission, which would have 
been filtered out by the background subtraction, or weak com- 
pact emission with sizes smaller than that of the APEX beam. 
Since methanol masers are almost exclusively associated with 
high-mass star forming regions it is very likely that the majority 
of these masers are located at the far side of the Galaxy. These 
masers have peak fluxes between 0.64 and 3.99 Jy beam -1 , well 
above the MMB surveys sensitivity of 0.17 Jy beam -1 , which 
would suggest that the methanol maser surveys can reveal some 
high-mass star forming regions that are missed at the sensitivity 
of ATLASGAL. 

A detailed analysis of the correlation and anti-correlation be- 
tween the ATLASGAL and MMB catalogues and their derived 
properties is beyond th e scope of this paper, b ut will be presented 
in a subsequent paper ( lUrquhart et al.L [2012b . 

6.3. Ultra-compact Hn regions 

We used the unbias ed 5 GHz VLA radio continuum survey by 
iBecker et al. ( 1994) to study the association between UCHn re- 
gions and submm clumps. In this VLA survey, the Galactic range 
from I = 350 - 40° and \b\ < 0.4° was covered with the C and 
CnD configurations, resulting in an angular resolution of ~4". 
This radio catalogue contains ~1300 sources, 30% of which they 
tentatively identify as UCHn regions. 

Using a search radius of 20" we find in the overlap region 
of both surveys 170 compact ATLASGAL sources that are posi- 
tionally coincident with a radio continuum emission (~6% of the 
ATLASGAL sample). Turning to the radio sources we find that 
23% of these are associated ATLASGAL sources and are there- 
fore likely to b e genuine UCHn reg ions, which is in reasonable 
agreement with lBecker et al.l (1 19941) initial estimate. 

6.4. Infrared Dark Clouds 

A considerable fraction of submillimetre br ight sources ap- 
pear d ark in Spitzer infrared images. Using the iPeretto & Fuller] 
(120091) catalogue of infrared dark clouds (IRDC) that is based 
on the Spitzer GLIMPSE images, we determined the association 
of ATLASGAL clumps with IRDCs. An association was based 
on overlap of the source ellipses in the two catalogues. In the 
part of the Galaxy common to both catalogues, we find that 30% 
of the IRDCs are associated with ATL ASGAL clump s . A sim i- 
lar analysis was recently conducted by I Wilcock et al.l (1201 2t) in 
a complementary part of the Galaxy {( = 300-330°) who deter- 
mined how many of the IRDCs were seen in the Herschel Hi- 
GAL 250/500 //m bands. They found a fraction of 38% which is, 
given the higher mass sensitivity of the Herschel survey, consis- 
tent with the ATLASGAL result. 

The IRDC peak optical depths are plotted against the peak 
flux densities measured in ATLASGAL in Fig.[T3] Since both 
quantities are proportional to column density, a correlation is 
expected. The main uncertainty in the determination of col- 
umn densities from ATLASGAL data is the dependence of the 
flux on the dust temperature, which might vary by about a fac- 
tor 2 (based on ammonia temperatures of ATLASGAL sources, 
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Fig. 13. Comparison between peak submm flux densities and 
8/iin peak optical d epths for ATLASGAL sources associated 
with IRDCs from the lPeretto & Fuller! (120091) catalogue. 
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IWienen et al.ll2.Q12b . Large uncertainties in the determination of 
IR optical depths result from the estimation of the foreground 
emission and the subtraction of the background. In addition, for 
high optical depths, th e extinction in the IR will saturate (e.g. 
IVasvunina et all 120091) . Given these major uncertainties of the 
latter method, column density estimates by optically thin dust 
emission in the (sub)millimetre still provide the most reliable 
method to measure the high column densities towards high-mass 
star forming clumps. 



6.5. MSX 21 fim association 

To study the mid-infrared properties of the ATLASGAL com- 
pact sources in more detail, w e extracted t he 21 /mi flux den- 
sities from the MSX survey dPrice et all 1200 ll) towards the 
ATLASGAL peaks. For such an analysis, the MSX data offer 
two advantages compared to the more sensitive and higher an- 
gular resolution data from the Spitzer MIPSGAL 24 fim survey 
(Carey et al. 2009): (i) the spatial resolution of the MSX data 
(18.3") matches the ATLASGAL resolution extremely well and 
(ii) the MSX data do not suffer from saturation that MIPSGAL 
shows towards strong sources, so that the full range of MIR 
brightnesses can be included in the analysis. This comparison 
is also more complete than a comparison with the MSX point 
source catalogue alone, since many ATLASGAL sources are as- 
sociated with more complex MIR emission not covered in the 
MSX point source catalogue. 

In Fig. [14] we show the correlation of 870 fim peak flux den- 
sity of the ATLASGAL sources and the 21 fim MSX fluxes mea- 
sured at the same position. To remove the contribution from the 
variable background seen in the 21//m images we have aver- 
aged the emission in an annulus centred on the ATLASGAL 
position with inner and outer radius of 30" and 120", respec- 
tively, and subtracted this from the 21 fim flux measured towards 
the centre of the ATLASGAL sources. The MSX fluxes span 
a large range and only a loose correlation of increasing MSX 
with ATLASGAL flux is seen. In fact, a large fraction of the 
ATLASGAL sources (~50%) show no 21 fim emission above 
the local background level, confirming again that many of the 



ATLASGAL compact sources do not have embedded or associ- 
ated luminous infrared sources. However, the fact that approxi- 
mately half of the ATLASGAL sources have MSX 21 fim emis- 
sion would suggest that the lower limit for the number associated 
with star formation is closer to 50%. This fraction is significantly 
higher than found from the straight comparison with the IRAS 
and MSX PSCs preformed in Sect. 6.1. This fraction is likely to 
increase when the ATLASGAL catalogue is compared with the 
Spitzer GLIMPSE and MIPSGAL surveys, which will be pre- 
sented in a future paper. 

These background corrected MSX 21 fim fluxes can be used 
to compare the mid-infrared properties of different subsamples 
ATLASGAL sources. In general terms, the ATLASGAL flux is 
a measure of the total mass of the clumps while the 21 fim MSX 
flux probes a warm component (~100K) of their spectral energy 
distributions. Therefore the ratio of both fluxes is a crude mea- 
sure of the fraction of heated material. We show the histogram of 
this ratio in Fig.[l5]and compare the distribution of this ratio with 
the corresponding hist ograms for sources from the RMS survey 
dUrquhart et al.Ll2008l) . both for identified high-mass young stel- 
lar objects and HII regions, and with the histogram of this rati o 
for CH 3 OH masers from the MMB survey dCaswell et al.Ll2010t) . 
The peak of the distribution for the RMS YSOs and HII regions 
is very comparable (green and blue histograms, respectively), 
and it is clearly offset from the ATLASGAL peak, again show- 
ing that with ATLASGAL many clumps with either no star for- 
mation or the very early stages of star formation are probed. The 
21 fim/870fim flux ratio for the methanol masers has a broader 
distribution (red histogram) which is consistent with covering 
the flux ratios from both, the bulk of the ATLASGAL and the 
RMS sources. 



7. Summary 

The ATLASGAL has mapped 420 deg 2 of the inner Galaxy at 
870jum, which traces thermal emission from dense dust clumps. 
Here we present a compact source catalogue covering 153 deg 2 
between Galactic longitude 330° < t < 21° and \b\ <1.5°. The 
source extraction algorithm SExtractor has been used to iden- 
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Table 3. Summary of associations between ATLASGAL com- 
pact sources and infrared and submillimetre surveys. 



Catalogue 


Coverage within ATLASGAL 


N sources 


IRAS 


Full survey 


110 


MSX 


Full survey 


1669 


IRAS-MSX 


Full survey 


645 


MMB 


-15° < / < 20°;|fc| < 2° 


459 


UCHii regions 


10° < / < 21°;\b\ < 1° 


170 


MSX 21/jm* 


Full survey 


3410 



* Number of ATLASGAL sources associated with MSX 21 fim emis- 
sion above the local background emission (see Sect. 6.5 for details). 



tify clumps and determine their parameters. To mitigate the vary- 
ing r.m.s. noise across the fields we performed the source de- 
tection on background subtracted signal-to-noise map, however, 
the source parameters were determined from the original emis- 
sion maps. In total 6639 sources have been identified above a 3<x 
threshold (~200mJy beam -1 ) and we find the catalogue is 99% 
complete above 6cr, which corresponds to ~0.4 Jy beam -1 . 

Comparing this catalogue of dense clumps to the MSX and 
IRAS point source catalogues we find ~40% are associated with 
far- and mid-infrared compact sources. The infrared emission 
likely arises from embedded young stellar objects (YSOs) or 
UCHii regions and is a strong indication that star formation is 
currently underway and is in an evolved phase in these clumps. 
This is supported by the positional correlation of 5 GHz ra- 
dio continuum emission with approximately 6% of these mid- 
infrared associated sources. In addition to the mid-infrared and 
radio associations we find a further ~500 clumps are associated 
with methanol masers (~7%). Methanol masers are thought to 
trace an earlier evolutionary stage than the mid-infrared emis- 
sion and have been almost exclusively found to be associated 
with high-mass star formation regions. A summary with the 
number of sources with associations among the different surveys 
is shown in Table[3] 



Combined the infrared, radio continuum and methanol maser 
associated sources make up approximately 50% of the cata- 
logue. These include all of the important evolutionary stages 
from young deeply embedded protostars identified by methanol 
masers to mid-infrared bright YSOs and UCHii regions and thus 
allows for the first time these different samples to be placed into 
a global framework. However, possibly more intriguing is the 
other 50% of the sample towards which none of the star forma- 
tion tracers have been detected as these sources are likely to be in 
an early stage of evolution. These sources are excellent targets 
for follow up observation, e.g. in the near future at very high 
spatial resolution with ALMA. 
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